About the lifetime of a bouncing droplet 
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When a droplet is gently laid onto the surface of the same liquid, it stays at rest for a moment 
before coalescence. The coalescence can be delayed and sometimes inhibited by injecting fresh air 
under the droplet. This can happen when the surface of the bath oscillates vertically, in this case the 
droplet basically bounces on the interface. The lifetime of the droplet has been studied with respect 
to the amplitude and the frequency of the excitation. The lifetime decreases when the acceleration 
increases. The thickness of the air film between the droplet and the bath has been investigated 
using interference fringes obtained when the system is illuminated by low pressure sodium lamps. 
Moreover, both the shape evolution and the motion of the droplet center of mass have been recorded 
in order to evidence the phase offset between the deformation and the trajectory. A short lifetime 
is correlated to a small air film thickness and to a large phase offset between the maximum of 
deformation and the minimum of the vertical position of the centre of mass. 

PACS numbers: 47.55.D-, 68.03.Cd, 68.15.+e 



I. INTRODUCTION 

The manipulation of tiny amounts of liquid plays a 
key role in any micro-device used in microfluidic appli- 
cations. A droplet is often laid onto a substrate where it 
interacts with its environment. Some liquid may be lost 
onto the substrate or the separation with the surround- 
ing liquid may be difficult. The idea is to manipulate 
the droplet without touching it [ij or in other words, to 
transport the droplet using another fluid, for example air. 
Couder et al discovered that it is possible to make an oil 
droplet bouncing onto a vibrating bath of oil. Some air is 
squeezed between the droplet and the bath avoiding coa- 
lescence. The question about the stability of the bounc- 
ing droplet is naturally raised. The aim of the present 
work is to discuss the lifetime of the bouncing droplet by 
studying its trajectory and the dynamics of the air film 
located between it and the bath. 
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FIG. 1: (left) Experimental setup for the droplet trajectories 
and deformations characterizing composed by the electromag- 
netic shaker GW-V55, a plate on which the container filled 
with the oil bath is attached, an incandescence lamp and a 
high speed video camera, (rigth) In order to observe inter- 
ference fringes, the container is lit homogeneously by three 
sodium low vapor pressure lamps. The droplet is filmed from 
the top using a high speed video camera. 



When a droplet is gently laid onto the surface of the 
same liquid, it stays at the surface during the time that 
the air film located just below the droplet is drained out. 
As soon as the air film thickness reaches a critical value 
(about 100 nm), the Van der Waals forces become of im- 
portance and the air film eventually collapses. Some non- 
coalescing drop systems have also been studied. It has 
been shown that, when a liquid-vapour CO2 interface is 
shaken horizontally, droplets may be produced and set 
in motion on the interface 0- Sreenivas et al. trapped 
droplets on an hydraulic jump Couder etal. 0,1 i 
have discovered that when the interface is vertically vi- 
brated using an electromagnetic shaker (cf. Fig. [1]), the 
droplet may avoid the coalescence. Indeed, when the 
droplet bounces, the air film is renewed. The relevant 
parameter is the reduced acceleration F defined as the 
ratio between the maximum acceleration Auj'^ and the 
gravity constant acceleration t/, A and uj beeing respec- 
tively the amplitude and the pulsation of the oscillation. 
The bouncing can be obtained when F is larger than a 
critical acceleration Tn that scales as the square of the 
frequency / = uj/2'k [J|. This criterion allows to deter- 
mine whether the droplet bounces or not according to 
the couple of parameters A and w. 

That phenomenon is perturbated at high accelerations 
when Faraday waves appear. This high acceleration limit 
F p depends on the viscosity of the oil in the bath. Con- 
sidering a given droplet size and a given viscosity, the 
bouncing phenomenon may only occur in a delimited re- 
gion of the (A, /) diagram, i.e. in between Tc and Tp 
(see Fig. H]). 

A question remains about the stability of the bouncing 
droplet. Indeed, at a given frequency and for an accel- 
eration larger than Tc{ijj), the droplet may bounce for 
a time longer than a few days [4]. On the other hand, 
for a large frequency (100 Hz) and for F > Tc{i^), the 
droplet bounces only for a few seconds before coalescing. 
The aim of this work is to measure the lifetime of the 
bouncing droplets and to determine where, in the phase 
diagram {A^oj), they bounce the longest time. 
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Due to the large number of parameters : A, f, the 
radius of the droplet R, the viscosity of the oil r?o, some 
have been arbitrary fixed. The viscosity and the droplet 
size are fixed to rjQ = 50cSt and 2R = 1.5 mm diam- 
eter respectively. In these conditions, the area we ex- 
plored in the A-f diagram is represented in Fig. [2] where 
the investigated set of parameters are represented. For 
each point, the lifetime has been measured. Then the air 
film dynamics has been studied by an interference fringe 
analysis. Moreover, an optical tracking method has been 
used to determine the position and the deformation of 
the droplet with respect to the phase of the plate. 

After the description of the experimental setup (Sect. 
II), the lifetime analysis will be developped (Sect. III). 
Explanation will be found by studying the interference 
fringes (Sect. IV) as well as the trajectory and deforma- 
tion of the droplet (Sect. V). 



II. EXPERIMENTAL SETUP 

A container is filled with silicon oil (Dow Corning 200 
Fluid) having a viscosity of 50 cSt. The plate is shaken 
according to a vertical oscillation of amplitude A and of 
pulsation ui (Fig. [T]left). For each fixed frequency, the 
acceleration is tuned between Tc and Ff by changing the 
amplitude of the vibration. The droplets are produced by 
using a syringe and are gently dropped onto the surface of 
the bath. The lifetime is measured using a digital timer. 
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The position and the shape of the droplets have been 
studied using image analysis. The motion of the droplet 
is recorded from the side using a fast video camera 
(1000 fps). Typical snapshots of a droplet are shown in 
Fig. [3l The one on the right presents a droplet squeezed 
by the plate. During a bounce, the deformation of the 
lower part and the upper part of the droplet is asymmet- 
ric. The evolution of its upper half part is easily deter- 
mined by contour detection on images while the lower 
part of the droplet motion is not accessible part of the 
time. This difficulty occurs particularly when the droplet 
deforms the bath surface. However, as the total volume 
and the evolution of the upper part of the droplet are 
known, the evolution of its lower part may be deduced as 
soon as a hypothetical shape has been speculated. The 
lower part shape is modeled by a flat bottom surrounded 
by a quarter of a torus [7]. The shape is fitted to the 
experimental contour taking into account that the total 
volume must be conserved. An example of this shape 
(white countour) is illustrated on Fig. [31 On the other 
hand, when the droplet is entirely visible, its shape is 
oblong, see Fig. [3] (left). The lower part of the drop is 
supposed to behave like a half ellipsoid. Eventually, the 
trajectory of the centre of mass is deduced as well as the 
evolution of the aspect ratio Ar = h/L, where h and L 
are respectively the height and the width of the droplet 

(Fig. m- 

In order to investigate the air film dynamics, three 
low pressure sodium lamps have been placed around the 
droplet (Fig. [l]left). That monochromatic light allows to 
visualize the interference fringes produced by the thin air 
film located between the droplet and the bath. A high 
speed video camera placed vertically above the plate has 
been used to record the motion of these fringes. The 
number of fringes gives information about the thickness 
of the air film [ 




FIG. 2: Phase diagram amplitude-frequency for a droplet of 
1.5 mm of diameter and made up of by oil of 50cSt. Top 
and bottom curves correspond to the Faraday Tp and Couder 
Tc acceleration thresholds respectively (see text). The data 
points represent the experimental conditions that have been 
investigated in the present work in order to cover at best the 
region where the droplet can bounce. Two large data points 
correspond to the cases I and II described in the text. 



FIG. 3: Snapshots of a droplet : (left) entirely visible, (right) 
at its maximum of deformation. The experimental contour 
of the droplet is represented by a white curve. The height h 
and the width L are indicated by the double arrows while the 
centre of mass is shown by the white circle. 



3 



III. LIFETIME OF DROPLETS 

Although Couder's criterion (F > Tc) is respected, 
we observed that bouncing droplets do not bounce for 
ever. Droplets lifetimes r give us information on their 
stability beyond that threshold. Avereged lifetimes are 
represented in Fig. 3] with respect to the reduced ac- 
celeration F. The different symbols correspond to the 
considered frequency, i.e. the squares, circles, triangles 
and rhombuses are for 25, 50, 75 and 100 Hz respectively. 
The vertical dashed lines represent the Couder thresholds 
according to the four different considered frequencies. 
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FIG. 4: Avereged lifetimes r as a function of the reduced ac- 
celeration r. The symbols correspond to different considered 
frequencies (see legend). The vertical dashed lines represent 
the Couder thresholds according to the different frequencies. 
The continuous line is the fit from Eq. [1] The large points 
correspond to the case I and II described in the text. 

The lifetime is found to drastically decrease with the 
reduced acceleration. One should note that, at 25 Hz, 
the lifetime is far below its real value. Indeed, the ex- 
periments were stopped when r > 1000 s. The droplet is 
qualified of stable beyond this arbitrary limit. A power 
law has been fitted (continuous line in Fig. |4]), namely 



1 



(F-l)i-' 



(1) 



A divergence is expected at F = 1, suggesting stability 
at very low frequencies. 

As far as the results corresponding to a given frequency 
are concerned, the maximum lifetime is found at the 
Couder threshold. Finally, the lifetimes corresponding to 
a same set of parameters are distributed along a broad 
distribution. This suggests that the mechanism explain- 
ing the life's duration of the droplet is probabilistic. 



In the next sections, two factors will be considered in 
order to explain the instability of the bouncing droplets. 



IV. INTERFERENCE FRINGES 




FIG. 5: Interference fringes obtained when the air film located 
between the droplet and the bath is the most squeezed; the 
number of fringes is maximum. The different pictures have 
been obtained for the four frequencies (a) 26 Hz, (b) 46 Hz, 
(c) 74 Hz, and (d) 98 Hz. 

One would expect that the instability of a bouncing 
droplet is due to the speed of the air drainage located 
between the droplet and the bath surface. In order to 
test this hypothesis, the bouncing droplet has been lit by 
sodium low pressure lamps. As this light may be consid- 
ered monochromatic, interference fringes can be observed 
through the air gap separating the droplet and the bath. 

The interference fringes produced by the air film have 
been recorded with a high speed video camera. During a 
certain period, the direct observation of the fringes shows 
that the circular fringes are visible when the droplet de- 
forms the surface of the bath. The fringes get closer 
and closer before separating and disappearing when the 
droplet takes off. When a movie taken at the beginning 
of the bouncing is compared to another taken at the end 
of the bouncing life, no evolution is observed. The mo- 
tion of the fringes remains periodical and identical dur- 
ing the whole life of the bouncing droplet whatever the 
frequency ! That observation means that the air film 
minimum thickness does not evolve. The death of the 
droplet cannot be explained by a progressive drainage of 
the air film. Taking this periodicity as an advantage, the 
observation frequencies have been chosen to obtain the 
maximum resolution time. Indeed, since the motion is 
periodical, an aliasing technique can be applied. 
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The movies allow to find the moment when the air 
film is the most squeezed, i.e. the phase ipc at which the 
fringes are the closest. In Fig. [5l four pictures corre- 
sponding to this instant are compared for f — 26 Hz, (b) 
46 Hz, (c) 74 Hz, and (d) 98 Hz. In Fig. 6, the phase 
(j)c is represented versus the acceleration F. The phase 
is calculated considering that the plate position follows 
a sinusoidal law Asin{u!t). Two significant observations 
can be made, (i) The number of fringes decreases when 
the frequency increases. Their number is related to the 
thickness of the film {3] , a small number of fringes is char- 
acteristic of a thin film. That means that the minimum 
thickness of the air film is smaller at high frequency than 
at low frequency, (ii) The phase (j)c increases with the 
frequency and the acceleration (see Fig. [6|) and seems to 
saturate towards 7r/2, i.e. when the plate reaches its max- 
imum position. In the next section, an additional study 
will be performed between the maximum of squeezing 
and the trajectory of the bouncing droplet. 
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FIG. 6: Phases <j!>c when the fringes are the closest as a func- 
tion of the reduced acceleration F. The different frequencies 
considered have been represented by different symbols (see 
legend). The dashed curve is a guide for the eyes. 



The probability of an air film to collapse is supposed 
to depend on its thickness. An activated process like 
an Arrhenius law may be imagined [l^, [HI . Under this 
hypothesis, the lifetime is determined by the minimum 
thickness that the air film reaches during one period. 
In order to explain why the minimum air film thickness 
changes with the frequency and the acceleration, the tra- 
jectory of a bouncing droplet will be analyzed in the next 
section. 



V. TRAJECTORIES AND DEFORMATIONS 

In this section, the droplet deformations are finely 
studied for two very different cases : (i) case I : a bounc- 
ing droplet characterized by a lifetime of about 420 s 
{A = 200 /im and / = 50 Hz) and (ii) case II : a bounc- 
ing droplet characterized by a lifetime of about 10 s 
{A = 250 /im and / — 100 Hz). Both cases are repre- 
sented by large symbols in Fig. 2 and Fig. 4. 




FIG. 7: Comparison of (a) the dimensionless position z/A of 
the centre of mass and (b) the aspect ratio Ar variation with 
respect to the reduced time (</? = uit) for cases I and II (see 
text). Cases I and II are represented by • and ♦ respectively. 
The continuous curve gives an indication concerning the phase 
of the plate. 

In Fig. [71 cases I (•) and II (♦) are compared along 
two parameters, (i) Fig. [7^ compares the positions of 
the centre of mass while (ii) Fig. [Th compares the varia- 
tions of the aspect ratios Ar with respect to the phase ip 
defined hy if = ujt where t is the time. Let us remember 
that a value ^4^ = 1 corresponds to a spherical droplet. 
The vertical position z of the centre of mass has been 
normalized by the amplitude A of the excitation. The 
continuous curves represent the position of the plate in 
arbitrary units in order to visualize its phase. The nor- 
malization of the time is justified by the fact that the 
motion is periodic from the birth till the death of the 
bouncing droplet. As for the interference fringes, no de- 
viation from these curves z/A((p) and Ar{ip) is observed 
during their life. 

Two particular phases will be discussed : (i) the phase 
'fimin at which the vertical position of centre of mass is 
minimum (see arrows Fig. [Tja) and (ii) the phase ifmax at 
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which the deformation is maximum, i.e. Ar is minimum 
(see arrows Fig. Hb). The normahzed position of the 
droplet is characterized by a larger amplitude in the low 
frequency case than in the high frequency one. However, 
the phase ipmm is nearly the same in both cases. These 
particular phases are indicated by arrows in Fig. [7^. It is 
noticeable that the trajectory in case I is more symetrical 
than in case II. 

The variation of the aspect ratio has nearly the same 
amplitude in both cases while the phases (pmax ai'e how- 
ever very different. The bouncing droplet is deformed 
during a smaller period in case I than in case II. Indeed, 
the aspect ratio Ar in case I is close and even above the 
unity during a longer period of time. This means that 
the time of interaction between the plate and the droplet 
is smaller than in case II. In that later case, the droplet 
is continuously deformed by the plate. 

In view of these results, a long lifetime is characterized 
by a large amplitude motion of the center of mass and 
a small phase difference between the maximum deforma- 
tion and the minimum of the center of mass. This phase 
difference can be explained by considering the inertia and 
the viscosity of the droplet. When the forced deforma- 
tion rate is too large (for high frequency), the droplet 
cannot follow the motion of the plate. A delay exists be- 
tween the motion of the plate and the deformation that 
is slowed down by the inertia and the viscosity. It re- 
sults in a stronger squeezing of the air film as observed 
in fringe figures. The ideal situation (for a long lifetime) 
is obtained when ipmin ~ y'max- 

VI. CONCLUSIONS 

We have studied the stability of bouncing droplets 
made up of sihcon oil of 50 cSt. The size of the droplet 



is 1.5 mm. For this system, the lifetime of a bouncing 
droplet has been found to decrease with F and to diverge 
for F = 1 and / = 25 Hz. Moreover, the distribution of 
the lifetimes for the same set of parameters suggests a 
probabilistic mechanism for the coalescence moment of 
the bouncing droplet. 

The motion of the air film and of the droplet has been 
tracked and recorded using image analysis. Both motions 
have been found to be periodical and identical during 
the whole life of the bouncing droplet. Such periodical 
motions indicate that the speed of the drainage cannot 
explain the decrease of the lifetime at high frequencies. 
On the other hand, the interference fringe analysis al- 
lows to determine that minimum thickness of the air film 
located between the droplet and the bath is smaller for 
high frequencies than for smaller frequencies. The film 
is more fragile at high frequencies and has consequently 
a larger probability of collapsing each time that the air 
film is squeezed by the plate. 

As for the measurement of the trajectory, the phase 
difference between the maximum of deformation and the 
minimum of the position of the center of mass has been 
found to be a key explaining the squeezing difference at 
low and high frequencies and therefore the difference of 
lifetime between the droplets. When this dephasage is 
close to zero and when the amplitude of the motion is 
large, the lifetime is increased. At high frequency, the 
phase difference increases and the interaction between 
the plate and the droplet is increased; the air film thick- 
ness is then decreased. 
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